Localization of energy for a Kerr black hole 
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Abstract 

In the teleparallel equivalent of general relativity the energy density of asymptot- 
icaly flat gravitational fields can be naturally defined as a scalar density restricted 
to a three dimensional spacelike hypersurface S. The scalar density has a simple 
expression in terms of the trace of the torsion tensor. Integration over the whole 
X yields the standard ADM energy. Here we obtain the formal expression of the 
localized energy for a Kerr black hole. The expression of the energy inside a surface 
of constant radius can be explicitly calculated in the limit of small a, the specific 
angular momentum. Such expression turns out to be exactly the same as the one 
obtained by means of the method proposed recently by Brown and York. 



PACS numbers: 04.20. Cv, 04.20.Fy 
(*) e-mail: wadih@guarany.cpd.unb.br 



I. Introduction 

Although there is a strong belief that Einstein's equations describe the dynamics of 
the gravitational field, it has not been possible so far to arrive at a definite expression 
for the gravitational energy in the context of Einstein's general relativity. On the one 
hand there are prejudices based on the principle of equivalence, which is often invoked 
to assure that the gravitational energy cannot be localized |TJ, 0. This principle, however, 
is not widely accepted as an essential element of the mathematical structure of vacuum 
general relativity 0. On the other hand, attempts based on the Hilbert-Einstein action 
integral fail to yield an unambiguous expression for the gravitational energy. The lat- 
ter is normally associated with surface terms in the action or in the Hamiltonian, and 
these surface terms do not exhibit the appropriate transformation properties under co- 
ordinate transformations (see the clear discussion in || in the case of asymptoticaly flat 
gravitational fields), a fact which prevents the construction of localized energy density. 

Recently an expression for quasi- local energy has been proposed by Brown and York||. 
Such expression is derived directly from the action functional A c \. The latter is identified 
as Hamilton's principal function and, in similarity with the classical Hamilton- Jacobi 
equation, which expresses the energy of a classical solution as minus the time rate of 
the change of the action, the quasilocal gravitational energy is identified as minus the 
proper time rate of change of the Hilbert-Einstein action (with surface terms included). 
Expressions for the quasilocal energy have been obtained for the Schwarzschild solution 
and for the Kerr solution 0. 

Einstein's equations can also be obtained from the teleparallel equivalent of general 
relativity (TEGR). The latter ammounts to a formulation which is geometrically different 
from the standard one based on the Hilbert-Einstein action integral, but whose dynamical 
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content is the same. In the TEGR the metrical quantity is the tetrad field and the 
action principle does not require the addition of non-covariant surface terms in the case 
of asymptotically flat gravitational fields. Because of this property, a localized energy 
density e(x) can be naturally defined from the action principle of the TEGR as minus 
the variation of the action with respect to the proper time N(x). Integration of e(x) over 
the whole three-dimensional space yields the ADM energy. Moreover, e(x) also appears 
in the expression of the Hamiltonian constraint C of the TEGR, a fact which allows the 
integral form of C = to be written as H — Eadm — 00, ||. 

We have calculated the energy inside a sphere of radius r Q in a Schwarzschild space- 
time by means of The expression turns out to be exactly the same as the one 
obtained through the procedure of ref.|| (expression (6.14) of 0). In this paper we 
consider the Kerr black-hole. We obtain the formal expression for the total energy in 
terms of non-trivial integrals in the angular variable 6. In the limit of slow rotation (small 
specific angular momentum) the energy contained within a surface of constant radius r Q 
can be calculated. Again the result obtained here is exactly the same as that obtained 
by Martinez who adopted Brown and York's procedure. The advantage of our proce- 
dure rests on the fact that the localized energy associated with a Kerr spacetime can be 
calculated in the general case, without recourse to particular limits, at least by means 
of numerical integration, whereas in Brown and York's procedure one has to calculate 
the subtraction term e° and for this purpose it is necessary to embed an arbitrary two 
dimensional boundary surface of the Kerr space S in the appropriate reference space (E 3 , 
say), which is not always possible 0. 

In section II we present the mathematical preliminaries of the TEGR, its Hamiltonian 
formulation and the expression of the energy for an arbitrary asymptoticaly flat spacetime. 
In section III we carry out the construction of triads for a three dimensional spacelike 
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hypersurface of the Kerr type, obtain the general expression of the energy contained in 
a volume V of space and provide the exact expression of the latter in the limit of slow 
rotation. Comments and conclusions are presented on section IV. 

Notation: spacetime indices ... and local Lorentz indices a, b, ... run from to 3. 
In the 3+1 decomposition latin indices from the middle of the alphabet indicate space 
indices according to [i — 0, i, a = (0),(i). The tetrad field e% and the spin con- 
nection uj^ab yield the usual definitions of the torsion and curvature tensors: R a = 
d^uj u a b + c uj u c b — T a = d^e a „ + a b e b „ — .... The flat spacetime metric is 
fixed by r/( )(o) = -1- 



II. The TEGR in Hamiltonian form 

In the TEGR the tetrad field e a M and the spin connection u>^ ab are completely inde- 
pendent field variables. The latter is enforced to satisfy the condition of zero curvature. 
The Lagrangian density in empty spacetime is given by 

L(e,co,X) = -ke(-T abc T abc + ^T abc T bac - T a T a ) + e\ ab ^ R ab , u (uj) . (1) 

where k = j^q, G is the gravitational constant; e = det(e a \ ab ^ u are Lagrange 
multipliers and T a is the trace of the torsion tensor defined by T a = T b ba - 

The equivalence of the TEGR with Einstein's general relativity is based on the identity 

eR(e,cu) = eR(e) + e(^T abc T abc + T abc T acb - T a T a ) - 2d,{eT») , (2) 

which is obtained by just substituting the arbitrary spin connection cu^b = °uj^ a b{e) + 
K^ a b in the scalar curvature tensor R(e,u) in the left hand side; u) lia b{e) is the Levi- 
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Civita connection and K^ ab — ~e a x ei, u (Tx fiV + T uXfl — T^ uX ) is the contorsion tensor. The 
vanishing of i? a ^(u;), which is one of the field equations derived from (1), implies the 
equivalence of the scalar curvature R(e), constructed out of e% only, and the quadratic 
combination of the torsion tensor. It also ensures that the field equation arising from the 
variation of L with respect to e a M is strictly equivalent to Einstein's equations in tetrad 
form (we refer the reader to refs.[^, § for additional details). 

It is important to note that for asymptoticaly flat spacetimes the total divergence 
in (2) does not contribute to the action integral. Therefore the latter does not require 
additional surface terms, as it is invariant under coordinate transformations that preserve 
the asymptotic structure of the field quantities^]. In what follows we will be interested 
in asymptoticaly flat spacetimes. 

The Hamiltonian formulation of the TEGR can be successfully implemented if we 
fix the gauge u>o a b = from the outset, since in this case the constraints (to be shown 
below) constitute a first class set||. The condition co> a6 = is achieved by breaking 
the local Lorentz symmetry of (1). We still make use of the residual time independent 
gauge symmetry to fix the usual time gauge condition e^) = e( )j = 0. Because of 
^Oab = 0, H does not depend on p kab y the momentum canonically conjugated to coW>- 
Therefore arbitrary variations of L = pq — H with respect to P kab yields Ukab = 0. Thus 
in view of u 0a b = 0, Ukab drops out from our considerations. The above gauge fixing can 
be understood as the fixation of a global reference frame. 

Under the above gauge fixing the canonical action integral obtained from (1) becomes 




(3) 
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H = NC + N*C t + S mn n m " + J-d k (NeT k ) + d k {W k Nj) . (4) 

07TG 

iV and N l are the lapse and shift functions, and S mri = — E nm are Lagrange multipliers. 
The constraints are defined by 



C = dj(2keT j ) - keZ kij T kij - ^(n%* - ^IT 2 ) , (5) 



C k = -e^d^ - n^% )ik , (6) 

with e = det(eQ) k ) and T l = g lk e^ l T^y k . We remark that (3) and (4) are invariant under 
global SO (3) and general coordinate transformations. 

We assume the asymptotic behaviour e^) k ~ r}j k + \hj k {^) for r — > oo. In view of the 
relation 

^ | tPxd^eTt) = ^ J s dS k (d t h tk -d k h tl ) = E ADM (7) 

where the surface integral is evaluated for r — > oo, we note that the integral form of the 
Hamiltonian constraint C = may be rewritten as 

J d?x\ke^T kl3 + ^(n%* - ^n 2 ) j = E ADM . (8) 

The integration is over the whole three dimensional space. Given that dj{eT^) is a scalar 
density, from (7) and (8) we define the gravitational energy density enclosed by a volume 
V of the space as 
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In similarity with Brown and York's procedure, we can also define the energy density as 
minus the variation of the action A?l with respect to the proper time N(x). Thus for a 
given set of solutions of the classical equations of motion the energy density e(x) can be 
defined as 

£ « = -WW) = • (10 > 

in agreement with (9). 



III. Energy of the Kerr geometry 

The Kerr solution Q describes the field of a rotating black hole. In terms of Boyer and 
Lindquist coordinates [|l(J (t, r, 9, </>) it is described by the metric 

ds 2 = -^[dt~asin 2 ed(t)\ 2 + ^f^[(r 2 + a 2 )d(j)-adt] 2 + ^-dr 2 + p 2 d9 2 , (11) 
p 2 P A 

A = r 2 — 2mr + a 2 , 

2 2,2 2n 

p = r + a cos ; 

a is the specifc angular momentum defined by a = — . The components of the metric 
restricted to the three dimensional spacelike hypersurface are given by gu = #22 = p 2 
and (733 = ^-sin 2 6, where £ is defined by 

S 2 = (r 2 + a 2 ) 2 - Aa 2 sin 2 9. 
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We define the triads as 



I -^sin6 coscf) pcosO cos<p —^sin9 sin(f>\ 

e {k)i = -^sinO sincf) pcosd sincj) ^sin8 costfi 

\ -j=cos(f) —psinO / 

(k) is the line index and i is the column index. The one form is defined by 



(12) 



e ( fc ) = e (fc) r dr + e d6 + e (fc) ^ , 

from what follows 

e« e(fe) = ^dr 2 + p 2 d9 2 + % S in 2 Qd<? 
A ' p 2 

We also obtain e = det(e^)i) = ~^ s ^ n ^- Therefore the triads given by (12) describe the 
components of the Kerr solution restricted to the three dimensional spacelike hypersurface. 

One readily notices that there is another set of triads that yields the Kerr solution, 
namely, the set which is diagonal and whose entries are given by the square roots of <fe. 
This set is not appropriate for our purposes, and the reason can be understood even in the 
simple clase of flat spacetime. In the limit when both a and m go to zero (12) describes 
flat space: the curvature tensor and the torsion tensor vanish in this case. However, for 
the diagonal set of triads (again requiring a — > and m — > 0), 



e M = dr , e ie) =rd6 , e w = r sinQ d<j> , 

some components of the torsion tensor do not vanish, T^)i2 = 1, ^(3)13 = sin6, and E g 
calculated out of the diagonal set above diverges when integrated over the whole space. 
Therefore the use of (12) is mandatory in the present context. 
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The components of the torsion tensor can be calculated in a straightforward way from 
(12). Only T( 3 ) 13 and T( 3 ) 2 3 are vanishing. The others are given by: 



o 

T(\\\2 = cosOcoscb (- H 1= sin 2 9 ==) 

K ' p pVA VA 



T (1)13 = sin9siri(f){ [2r(r 2 + a 2 ) - a 2 sin 2 9 (r-m)} + — + -^=} 

pS p 6 VA 



S AS 

7)i)23 = cosBsincp {p h a 2 sin 2 6 (— -)} 

p pS p A 



2 

V (X p 

T(2M2 = cosOsind ( — I == sin 2 6 7 =) 



1 7"5j 

n (2)i3 = —sinOcoscj) { [2r(r 2 + a 2 ) — a 2 sin 2 9 (r — m)] H H =] 

pS p 3 VA J 



S AS 

^(2)23 = -cosOcoscp {p ha 2 sin 2 6 1 (— -)} 

p pS p 4 



T (3 )i2 = sm#[-- + -^= + ^ 7 =cos 2 #] 
1 ' p \fK pVA 

In order to evaluate (9) we need to obtain T\ After a long calculation we arrive at 



! VA VA A 2 2 2-2/1/ M 

T = — - H ^r[2r(r + a ) - a sm 6>(r - m)\ , 

p z S p z S z 



m , . „ „ a 2 1 cosO S o . 2 A S, 

T = sinQcosQ— + — p h a 2 sin 2 6(— - 

p 4 pSsin# [F p v pS p 3 ' 



T d = . 
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The gravitational energy density inside a volume V of a three dimensional spacelike 
hypersurface of the Kerr solution can now be easily calculated. It is given by 



E n 



57T JV 



dr d9< 



d_ 

dr 



sin6[p 



E 
P 



p£ 



y 2r(r 2 + a 2 ) — a 2 sin 2 9{r — m) 



+ 



d_ 
89 



£a 2 . 2 _ _ cos 
sin Ucost) + 



Next we specialize E g to the case when the volume V is contained within a constant 
radius surface r = r Q , assuming r a to be greater that the outer horizon r + = m + 
\/m 2 — a 2 . The integrations in and r are trivial. Also, because we integrate 9 between 
and 7r, the second line of the expression above vanishes. We then obtain 



- [ d6 sinoi p H ( 2r(r 2 + a 2 ) — a 2 sin 2 9(r — m) 

4 Jo V p pS 1 



(14) 



We have not managed to evaluate exactly the integral above. However, in the limit 
of slow rotation, namely, when y- « 1 all integrals have a simple structure and we can 
obtain the approximate expression of E g . It reads 



E n 



To 1- 



\ 



2m a 2 \ » 

1 + ^7 + 



6r G 



2 + 



2m / 2m \ 
+ 1 + — 



2m a 2 " 1 



(15) 



This is exactly the expression found by Martinez pi for the energy inside the surface of 
constant radius r Q in a spacelike hypersurface of a Kerr black hole, in the limit of small 
specific angular momentum. As in ref. 0, we have not expanded the square root which 
appears in (15) in powers in 

Let us mention finally that the expansion of p + — in the integrand of (14) yields 
— So, whereas the remaining term corresponds exactly to e, expressions (3.17) and (3.1) 



respectively of ||. It does not seem to be possible, however, to split <9j(eT J ) into two 
terms such that their integrals arise in the form e — e . 

Comments 

The gravitational energy E g defined by (14) can be evaluated for an arbitrary value 
of a by means of numerical integration. This is the major advantage of our procedure 
as compared to that of Brown and Yorkj5|]. By means of the latter one cannot construct 
expressions like (13) and (14), which may be useful in the study of astrophysical prob- 
lems, since in a general situation Brown and York's procedure requires the embedding of 
an arbitrary two dimensional boundary surface of the Kerr space in the reference space 
E 3 , a construction which is not possible in general (the evaluation of Eq in || is only 
possible in the limit f- « 1). Therefore the present approach is more general than that 
of ref. ||. Finally we remark that we expect expressions (9) and (10) to be useful in the 
study of the thermodynamics of self-gravitating systems, where the gravitational energy 
plays the role of the thermodynamical internal energy that is conjugate to the inverse 
temperature. We hope to come to this issue in the future. 
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